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ABSTRACT 

Using HST NICMOS 2 observations we have measured 1.6 /um near infrared nuclear luminosities 
of 100 3CR radio galaxies with z < 0.3, by modeling and subtracting the extended emission from 
the host galaxy. We performed a multi-wavelength statistical analysis (including optical and radio 
data) of the properties of the nuclei following classification of the objects into FR I and FR II, and 
LIG (low-ionization galaxies), HIG (high- ionization galaxies) and BLO (broad-lined objects) using 
the radio morphology and optical spectra, respectively. The correlations among near infrared, optical, 
and radio nuclear luminosity support the idea that the near infrared nuclear emission of FR Is has 
a non-thermal origin. Despite the difference in radio morphology, the multi-wavelength properties of 
FR II LIG nuclei are statistically indistinguishable from those of FR Is, an indication of a common 
structure of the central engine. All BLOs show an unresolved near infrared nucleus and a large near 
infrared excess with respect to FR II LIGs and FR Is of equal radio core luminosity. This requires the 
presence of an additional (and dominant) component other than the non-thermal light. Considering 
the shape of their spectral energy distribution, we ascribe the origin of their near infrared light to hot 
circumnuclear dust. A near infrared excess is also found in HIGs, but their nuclei are substantially 
fainter than those of BLO. This result indicates that substantial obscuration along the line-of-sight to 
the nuclei is still present at 1.6 /im. Nonetheless, HIGs nuclei cannot simply be explained in terms of 
dust obscuration: a significant contribution from light reflected in a circumnuclear scattering region 
is needed to account for their multiwavelength properties. 

Subject headings: Galaxies: active - Galaxies: elliptical and lenticular, cD - Galaxies: nuclei - Galax- 
ies: evolution - Galaxies: photometry - Galaxies: structure - Infrared: galaxies 



1. INTRODUCTION 

Infrared observations of radio galaxies are a useful tool 
to explore the physics of their nuclei in the framework 
of the Unified Models for Active Galactic Nuclei (AGN) 
as in this band the impact of dust absorption is strongly 
reduced with respect to the optical. In particular, the 
study of the infrared nuclear sources allows us to inves- 
tigate the 'scheme' unifying different classes of AGN, by 
exploring the properties of their accretion disks and the 
presence of any absorbing material that can account for 
the complex AGN taxonomy. 

In the "zeroth-order approximation" of the AGN 
unification scheme for radio-loud sources (e.g. 
lUrrv fe Padovanil [1995), powerful ra dio galaxies with 
FR II edgc-brightcned morphology (Fanaroff & Riley 
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fl97l are believed to be misaligned quasars, while lower 
power, edge-darkened FR Is are associated with BL Lac 
objects. However, it is clear that this zeroth order 
approximation picture based on statistical comparison 
of the properties of radio-loud AGN is over-simplified. 
Inconsistencies between the characteristics of the 
"parent" and "beamed" populations are well reported 
in the literature: for exa mple, several BL Lacs show 
broad emission lines (e.g. IVermeulen et al.l 11995 ) ; and 



there are inconsis tent environments (e.g. IZirbel 119971 : 



lOwen etldl Il996j) and inconsist ent radio morphologies 
(e.g. lAntonuccilll986l) (see also lUrrv fc Padovanil 1 19951 
and IChiabergdl2004l for reviews on this subject). 

Furthermore, how the large-scale radio properties re- 
late to different nuclear properties is still not clear. 
In fact, there is mixed evidence for a correlation be- 
tween the morphological FR I/FR II dichotomy and 
different levels of nuclear activity (iBaum et all 119951: 
I Chiaberge et all 119991: IXu et al.lll999al IChiaberge et al.l 



2002a; IVerdoes Kleiin et al.l |2002j). Another clas- 



sification, b ased on the optical narrow emission 
lines ratios dHine fc Longairl fl979l lLaing et al I 11994 
Uackson fc RawlingsHl997t iButtiglione et alJ l2010) intro- 
duces two groups of radio-loud AGN: the Low Exci- 
tation Galaxies (LEG) and High Excitation Galaxies 
(HEG) 7 . Radio galaxies with FR I radio morphology 
almost always belong to LEG class, while nearly all 

7 In the following we will refer to those objects as High Ionization 
Galaxies (HIG) and Low Ionization Galaxies (LIG). See Sect. 2 for 
details. 
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HEGs are FR II radio sources. However, there is a 
group of FR II LEGs whose place in t he unifying sce- 
nario is sti ll largely not understo od. I Wall fe Jackson! 
(11997ft and IJackson fc Wall (j!999ft proposed that such 
objects constitute a single population of radio galax- 
ies together with FR Is. In line with this idea, the 
two spectroscopic classes, LEG and HEG, show signif- 
icant (fiff<erena3s_Jn^ formation 
(e.g. IBaldi fc Capettil 120081: lSmolcicll2009ft. env ironment 
(e.g. IChiaberge et alj|2000al lHardcastldl2004j ). trigger- 
ing mechanis m of nuclear activity (g as-rich or gas-poor 
merger) (e.g. IBaldi fc Capettil l2008l i: there are indica- 
tions that they differ also in terms of accretion rate 
and r adiative efficiency of dis k (e.g. iMarchesini et al.l 
120041: iBalmaverde et al.l |2008|) ; contrarily to what is 
seen in HEG, LEG do not show evidence for large 
absorbing column den sities in the X-ray spectra (e.g. 
lHardcastle et al.ll2006l). for warm dusty tori in mid-IR 
images ( |van der Wolk et "al"1l2009ft , and for broad lines in 
their optical spectra (e.g. iButtiglione et alll2010ft . Taken 
together these differences suggest that high and low ex- 
citation galaxies belong to different classes of radio-loud 
nuclear activity which are associated with distinct en- 
vironmental conditions. Furthermore, it has been sug- 
gested that in the framework of the unification model, 
HEG and LEG might represent the p arent population 
of QS Os and BL Lacs respectively (e.g. IJackson fc Walll 
fl999l) . 

A number of recent papers in the literature address 
the topic of t he infrared view of the active nucleus of ra- 
dio galaxies. iDicken et al.l (|2009[ ) show that the infrared 
emission at both 24 and 70 micron can be explained 
in terms of dust heated predominantly by the AGN. 
More comprehensive mid-IR and far-IR spectrophotom- 
etry for the z < 1 3CR objects has been obtain ed from 
ISO (jSiebenmorgen et a l. 2004; Haas et al. 2004) as well 
as from Spitzer fe.g.. IShi et~aLl 120051: lHaas et all 120051 : 
lOgle et al.|[200l IClearv et al.ll2007l: iLeipski et al.ll2009l ). 
For high-redshift 3CR sources near infrared spectropho- 
tomet ry us ing Spitzer has been presented by lHaas et al.l 
(2008) and lNetzer eTall (|200l . 

The mid-infrared observations are sensitive to thermal 
emission from warm dust and have revealed rather com- 
plex properties within the different radio-galaxy classifi- 
cations. The majority of FRIs are faint at those wave- 
lengths, implying that the bulk kinetic power of the ra- 
dio jets is the ene rgetically dominant mechanis m in most 
of these sources (jWhvsong fc Antonuccill2004ft. Even in 
FR Is this is not a universa l conclusion a s lHaas et aT] 
(2004) with ISO data and ILeipski et all (|2009l ) with 
Spitzer data show that in some FR Is a dust emis- 
sion b ump is detected. On the other hand, lOgle et all 
(2006) found that about 50% of the FR II narrow-line 
radio galaxies have mid-infrared photometric and spec- 
troscopic properties comparable to quasars of matched 
radio luminosity, while there is a subsample of mid-IR 
weak FR II radio galaxies which may constitute a sepa- 
rate population of non-ther mal, jet-dominated sources 
with low accretion power. lHaas et ail (|2004t ) present 
evidence that the MIR data do in fact support the 
orientation-dependent unification scheme of the power- 
ful FR II galaxies with quasars. This has been sub- 
sequen t ly confirmed by f urther work, e.g. | Haas et all 
(2008T): lold7et~aT1 pOfil ): IDicken et al.l (12009ft . 



This paper has the aim of investigating AGN activity 
by focusing on the nuclear 1.6 /im near infrared emis- 
sion, identifying its physical origin and relating it to 
both th e radio morphological and opt ical spectroscopic 
classes. IChiaberge et al.l ([19991 l2002aft studied the HST 
optical images of the nuclear regions of radio galaxies 
belonging to 3CR sample with the same purpose. The 
NIR snapshot survey of 3CR objects with z < 0.3 pro- 
vided us with a larger coverage of the sample than that 
in the optical band. Furthermore, the fraction of 3CR 
sources with spectroscopic classificat ion is now substan- 
tially larger (IButtiglione et al.ll2010ft . allowing us to ex- 
plore in greater detail the behavior of the nuclei of each 
class. 

The resolution of HST images is suitable for this pur- 
pose and allow us to distinguish the NIR emission of 
stellar host galaxy background from the genuine AGN 
radiation. 

The structure of the paper is as follows. In Sect. 2 we 
describe the sample and in Sect. 3 we derive the surface 
brightness profiles which we fit with either Sersic or core- 
Sersic models. The multiwavelength properties of our 
sample are discussed in Sect. 4 in order to identify the 
origin of the NIR nuclear emission tackled in Sect. 5. In 
Sect. 6 we summarize and draw conclusions. In Appendix 
we compare the FR I radio core flux observed with VLA 
and VLBI. 

We adopt a Hubble constant of Ho = 75 km s _1 
Mpc -1 and qo = 0.5. We used this cosmology in or- 
der to be in accordan c e with t he similar results o btaine d 
bv IChiaberge et all (11999ft . IChiaberge etail (|2000ah . 
IChiaberge et al.l (l2002al T Assuming a different cosmol- 
ogy would not affect our results since the redshift of the 
sources considered in this paper is limited to z = 0.3. 

2. SAMPLE AND HST OBSERVATIONS 

The sample of galaxies we analyzed b elongs to the Re- 
vised Third Cambridge Catalog 3 CR (jBennettl H962bl 
lBennettlfl962allSpinrad et al.lll985ft . Since it is selected 
based on the low-frequency radio flux (SirsMHz > 10 Jy), 
the sample is essentially free from orientation bias. We 
chose all 116 3CR sources with z < 0.3 with the primary 
goal of characterizing the radio galaxy hosts mostly free 
from the effects of dust. One hundred objects were ob- 
served as part of HST/NICMOS 3CR snapshot program 
G O 10173 (PI: Sparks) . The NIR images are p ublished 
in iMadrid et ail (120061 ) and iFlovd et alj (j2008h . Only 
eleven galaxies were observed as a part of other programs 
and details of these are given in TableHJ We also included 
in the sample the FR I radio galaxy NG C 6251 because, 
as pointed out bv lWaggett et al.l (|1977ft it should be in- 
cluded in 3CR catalog. 

All of the objects of our sample were observed with 
HST NICMOS Camera 2 (NIC2), which has a field of 
view of 19'.'2 x 19'.'2 and a projected pixel size of 0.075", 
using the F160W filter (similar to H band), which is cen- 
tered at 1.60 /im, covering a wavelength range from 1.4 
to 1.8 /im and which includes the Pa/3 line for objects 
with z > 0.1. Using the bro ad and narrow line Ha fluxes 
measured for 3CR galaxies (Buttigli one et al.ll2009ft over 
a 2" x 2" aperture, we have estimate the contamination 
from line emission in the H band. In the case of the nar- 
row line s we assumed rati o Pa/3/Ha = 0.06 from atomic 
physics (|Osterbrockl 11989ft . Since the broad line galax- 
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TABLE 1 
Log of observations 



Name 


Date 






yy-mm-dd 


S 


3C 084 


1998-03-16 


639.9 


3C 264 


1998-05-12 


447.7 


3C 270 


1998-04-23 


767.8 


3C 272.1 


1998-07-13 


767.8 


3C 274 


1997-11-20 


127.9 


3C 293 


1998-08-19 


2751.8 


3C 305 


1998-07-19 


2783.8 


3C 317 


1998-08-26 


639.8 


3C 338 


1997-12-18 


4159.8 


3C 405 


1997-12-16 


1343.9 


NGC 6251 


1998-07-06 


543.5 



Note. — Log of observations for the objects not included in the 
HST 1.6 lira near i n frared snapshot p rogram 10173 (PI: Sparks) 
UMadrid et al.l 120061 ; IFlovd et aU 120081) . All data are taken with 
the same configuration with Camera NIC2 and F160W filter. 



ies are know not to conform to case B we have used the 
observed near infra red broad line regi on rat ios observed 
in typ e 1 AGN by iRiffel et al.l (|2006[ ) and lLandt et al.l 
(2008) to estimate the contamination in the F160W fil- 
ter. This analysis shows the contamination to be < 15% 
for Broad Line Objects and < 4% for Narrow Line Ob- 
jects and in both cases that is small enough not to affect 
our conclusions regarding the origin of the nuclear NIR 
emission. 

Due to technical difficulties, 3C 410 and 3C 442 were 
only observed for 575.9 s, while all the remaining galaxies 
were observed for a total exposure time of 1151.8 s. A de- 
tailed descriptio n of the observations and data reduction 
can be found in lMadrid et~aT1 (|2006l ). 

Each galaxy was classified in two ways. Firstly, 
into FR I and FR I I based on radio morphology 
(|Fanaroff fc Rilevl Il974f ): the radio mor phological clas- 
sificat ion, mostly being adopted from IZirbel fe~ Baum 
ljl995t ). Secondly, based on the optical broad line 
strengths and narrow emis sion lines ratios (from 
iButtigl lone "et~aT1 [20091 l20inh we classified them into 
the stand groupings of High and Low Ionization Galax- 
ies (HIG and LIG) and Broad Line Objects (BLO). 
Note that the BLO objects belong to HIG class be- 
cause they show high-ionization narrow emission lines. 
We adopt the HIG/LIG nomenclature to better repre- 
sent the physical condition of the narrow line region 
gas in these objects. This classification is however 
entirely consistent with the HE G/LEG scheme widely 
adopted in the l it erature (e.g.. iHine fc Longairl 119791 : 
Laing et all Il994t iButtiglione et al.l \201(£ . Following 
Buttiglione et al.l ()2010f ) we introduced a 4th class of 
Extremely-Low Excitation Galaxies (ELEG), character- 
ized by an extremely low value of [O IIIJ/H/3 ~ 0.5, ^6 
times lower than usually observed in 3CR/LIG. In this 
work we consider ELEGs as a separate class of objects. 
The objects for which we d o not have new spect r al dat a 
are classified according to Uackson fc Rawlingsl (|1997l ). 
In the sample there is also one galaxy, 3C 198, charac- 
terized by a star forming spectrum (marked as "SF" in 
Table[2]) and a BL Lac object (3C 371) which are not con- 
sidered in the following statistical evaluations because of 
their unique properties among the objects of our sample. 
For one galaxy (3C 410) we could not find any reliable 
spectral classification in the literature. 



3. RADIAL BRIGHTNESS PROFILE FITTING AND 
NUCLEAR LUMINOSITY MEASURING 

The aim of this paper is to measure the NIR nu- 
clear sources. In order to disentangle such a compo- 
nent from the underlying stellar emission of the host 
galaxy, it is necessary to model the distribution of 
the galactic light. To this purpose we derived the 
surface brightness profile of all our sources. Radial 
profiles can be derived using a number of different 
methods: ellipse fitting to the isoph otes of the galax- 
ies (e.g., iHeraudeau fc Simienl I1996D obtainin g a one- 
dimensional analysis (e.g.. iBaggett et~aT1 H998h. an d full 
two-dimensional analysis (e.g.. iBvun fc Freeman! 119951 : 
iPeng et ail 120021 ). In this paper, we adopted the first 
approach. 

Th e "ellipse" task in IRAF STSDAS (jJedrzeiewskil 
Il987h is suitable to analyzing the axisymmetric structure, 
evident in most objects. This task allows us to trace the 
radial profile and radial changes in ellipticity, position 
angle, and coordinates of the center of the isophotes. 

Most of the host galaxies of 3CR sources appear as 
smooth ellipticals in infrared. Only in few cases are dis- 
tortions due to dust lanes, jets, or plumes. Before deriv- 
ing the galaxy profile, we mask out all spurious objects 
such as stars, jets, dust lanes, bright, galaxy companions 
and residual diffraction patterns associated with bright 
stars. When dust features are present in the nuclear re- 
gion we adopted a multi-stage strategy. To stabilize the 
fit in the radial range affected by the dust we started by 
fixing some or all fit-parameters to the values of the last 
isophote that is not affected by the dust. We then pro- 
duce an initial model image and we subtracted it from 
the original data to unveil the presence of further faint 
sources and/or dust features, that were then also masked 
out and so on until no further feature were found. 

In Fig. Q]we show the stages of the procedure to derive 
and to fit the radial brightness profile. 

3.1. Profile fitting 

The surface brightness profile can be fitted using a va- 
riety of functions. We used two different models, Sersic 
and core-Sersic. The former (|Sersidll968f ) is a law of the 
form: 



I(r) 



-br,[(^f -1] 



where j3 is the concentration parameter and /3 = 1/n, 
the inverse of the Sersic index n. I e is the intensity at the 
effective radius r e . For a Se rsic model with 5 < n < 10 
b ra » 1.9992n- 0.3271 (e.g.. ICapaccioh1 ll989l ICaon est al.l 
fl99ilGraham fc Driverll2005l) . 

The latter is composed of an outer Sersic model and 
a inner power-law profile, suitable to describe the in- 
nermost part of the profile when the object shows a 
deficit toward the ce nter with respect to the Sersic law. 
iTruiillo et all (|2004l ) argued that it is convenient to re- 
strict the core-Sersic model to an infinitely sharp tran- 
sition between the Sersic and power-law profiles at the 
break radius it, assuming this form: 



J(r) = 1(b) 



(7) *(rb-r) + e 
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Fig. 1. — This figure summarizes the fitting procedure of the surface brightness profile (in this example 3C 349). In the upper left 
panel we show the original NICMOS image. The lower right panel presents the model image using the parameters obtained by the task 
ellipse. The lower left panel shows the residuals image (typically ~10% in the nuclear region). The object in the S-W corner in the residual 
image was masked and not fitted by the ellipse task. All the images are ~7" wide. The upper right panel shows the fitting of the radial 



unidimensional brightness profile (in units of erg s 



arcsec ) with a Sersic law (blue dashed curve), a nuclear component 



(black dot-dashed line) and a background level (green horizontal line). We take into account in the fit the convolution (black solid line) of 
the model with the instrumental PSF. The final model is the sum of each component and is the red solid line. The vertical dashed line 
gives the location of the effective radius r e . In the lower region of the panel we show the residuals. 



where r;, separates the inner power law of logarithmic 
slope 7 from the outer part of the profile described by a 
Sersic model. 6*(r&-r) is the Heaviside step function. 

We fit the brightness profiles with both models by min- 
imizing the residuals, leaving free to vary the 5 (or 7) 
variables: the background level, the three parameters 
describing the Sersic (or five for core-Sersic) model and 
the nuclear component. Fig. [2] illustrates four examples 
of fits to the NIR brightness profiles. In Sect. I3.2l we de- 
scribe the method used to select the appropriate law for 
each galaxy. 

To take into account the PSF effects, we convolved 
the models prior to t he co mparison with the data 
iCapetti fc Balmaverdd (12001 . The PSF of NICMOS- 
NIC2 at the center of the chip was derived using the 
TINYTIM program written by J. Krist. The PSF image 
was then convolved with a model galaxy with a power 
law profile, which represents the innermost region of a 



Sersic/core-Sersic model. This was done for a set of 
galaxy models with different profile slopes 7', all of them 
having zero ellipticity. We then derived the radial profile 
of the convolved model galaxies using circular isophotes, 
and we calculated the intensity ratio between the con- 
volved model galaxy and the original model galaxy for 
different radii. With this method we obtained a set of 
radial ratios for different values of the innermost slope 
that can be applied to the models to "simulate" the two- 
dimensional convolution. 

The appropriate value of 7' was measured as part of 
the fitting procedure. In the case of a Sersic model, 7' 
is the value of the logarithmic slope of the profile at 
O'.'l, which approximately corresponds to the resolution 
of NICMOS-NIC2. In the core-Sersic model, 7' is the 
lo garithmic slop e of th e inner power-law profile. 
Donzelli et al.l (|2007| ) already performed the fitting of 



the radial brightness profiles of 3CR sample with Sersic 
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Fig. 2. — This figure shows four explicative fitting methods used. In the upper parts we report the NICMOS images of the 4 radio-galaxies. 
In the lower part each panel represents the modeling of the surface brightness profile with a Sersic or a core-Sersic. The description of lines 
are as in Fig.[l] We mark the position of the effective radius with a black vertical dot-dashed line and, for a core-Sersic model, the break 
radius with a dashed vertical line. More in details we report four examples: in the upper left panel 3C 17 was fitted with a Sersic model; 
in the upper right panel 3C 123 with a core-Sersic model to obtain an upper-limit to its nuclear luminosity; in the lower left panel 3C 129 
with a Sersic model in the case of an weak nuclear source; and in the lower right panel 3C 296 with a core-Sersic with a weak nuclear point 
source. In this last panel the nuclear source plotted in black dashed line is amplified of a factor 10 for a better visualization. In the plots 
of 3C 129 and 3C 296 the green horizontal line of the background level is absent because it is below the range of intensity shown. 



law, focusing on their large scale behavior, thus not in- 
cluding t he central region of th e galaxies. For the objects 
in which iDonzelli et al.l (|2007f ) found a disk-like compo- 
nent in the galaxy profile, we exclude from the fit the 
regions dominated by the disk. The pa rameter values 
found are mostly consistent with those of IDonzelli et al.l 
(2007) for the objects in which the nuclear component 
is not dominant with respect to the galaxy stellar emis- 
sion. The resulting models produce residuals in the form 
of large scale fluctuations with a typical amplitude of < 
5%, indicative of the goodness of the fit. 

Besides the parameters describing the host galaxy pro- 
file, the fitting procedure also yields the flux of the nu- 
clear source. 



To assess the reliability of the detection of a nuclear 
source we adopt the following operative approach. We 
extract the radial brightness profile using the IRAF "rad- 
prof" task and measure the nuclear FWHM, setting the 
background level at the intensity measured just beyond 
the nuclear point source at a radius of ~0.3". If the 
FWHM is consistent with that of the PSF (~ 0.14 arc- 
sec), we took this as an indication of the presence of 
an unresolved source. Alternatively, if there was no ev- 
idence for a nuclear point source, we proceed to fit the 
profile with a core-Sersic model leaving all parameters 
free to vary, but fixing 7 = 0. This corresponds to the 
most conservative choice for 7 as it minimizes the galaxy 
contribution and maximizes the nuclear source. We con- 
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sider the nuclear flux thus derived as the upper limit to 
the intensity of the central source. 

3.2. Sersic or core-Sersic? 

The choice between Sersic and core-Sersic model for a 
given galaxy is related to the corresponding minimized 
X 2 value obtained after the fitting procedure. Since 
the core-Sersic model includes more free parameters, the 
best core-Sersic model usually will reproduce better the 
galaxy profile than the Sersic model. Therefore, to de- 
termine the goodness of the fit, we use the "reduced" \ 2 
(Xr) — X 2 1 i n ~ P)i where n is the number of fitted data 
points and p is the number of fitted parameters, which 
can be either 3 or 5 for Sersic or core-Sersic models, re- 
spectively) . We consider the best model that corresponds 
to the Xr closest to the unity. Qualitatively, the choice 
can be justified by the form of the residuals. In fact, when 
fitting a core-Sersic galaxy with a pure Sersic model, the 
residuals in the central regi on show the characte ristic S- 
shape, already discussed bv lTruiillo et al.1 ()2004f ). 

Nevertheless, since our objects are at distances of up 
to z = 0.3, we can misclassify a genuine core galaxy as 
a pure Sersic due to an insufficient physical resolution of 
the image. Th e typical core size of an ellipt i cal galaxy is 
~200 pc (e.g. ICapetti fc Balmaverdd l2005t iFaber et al.l 
Il997t iFerrarese et al.l |2006[ ) which corresponds to 0'.'14 
at z = 0.15. Considering the objects at z > 0.15 with de- 
tected nuclei (excluding the BLOs, for which the nucleus 
outshines the host galaxy emission in the central region) , 
for which there is ambiguity on the presence of a core, 
the typical discrepancy between the estimates of the nu- 
clear luminosity obtained adopting either a Sersic law or 
a core-Sersic model is ~ 30 %, and at most a factor of 
2. The sources potentially involved in this ambiguity are 
only ~15% of the sample and thus they do not alter the 
main results of this paper. 

The presence of a bright point source may also limit our 
ability to see a possible shallow core and to discriminate 
between Sersic and core-Sersic models. However, in those 
cases (e.g. in the case of the BLOs, Fig|3j), the difference 
between the nuclear fluxes obtained adopting the two 
different models is irrelevant, being smaller than ~ 2%, 
thus well below the typical photometric error. 

4. RESULTS 
4.1. Detection of nuclei 

In Table[5]we report the 1.6 NIR nuclear fluxes and 
luminosities for the sample measured with the procedure 
explained in Sect. [3l We also give the optical (R band), 
radio (nuclear at 5 GHz and large scale at 178 MHz) 
fluxes, and [O III] emission line luminosities for all the 
objects of the sample. 

Here we examine the detection rate of the NIR central 
sources (Table [3]) , considering each radio-morphological 
and optical-spectral class separately. Overall, we detect 
nuclei in 64%±5% of the entire 3CR sample. We de- 
tect nuclei in 81%±8% and 58%±6% of the FR Is and 
FR lis, respectively. Among the FR lis, all of the BLOs 
have an unresolved nuclear point source, while the de- 
tection rate of nuclei in LIGs and HIGs is 44%±10% and 
54%±9%, respectively. The only one BL Lac included 
in the sample has a bright nuclear source. The nucleus 
of 3C 410, the only spectroscopically unclassified FR II, 







3C390.3 ~ 





-1.5 -1.0 -0.5 0.0 0.5 1.0 

log(r [arcsec]) 

Fig. 3. — The surface brightness profile of 3C 390.3 fitted with 
a Sersic model (green solid trace) and a core-Sersic model (red 
dot-dashed trace). The lines below represent the two models used, 
which include the PSF-convolved galaxy model, the nuclear source 
and the background (the green solid and the red dot-dashed lines 
correspond to the Sersic and core-Sersic, respectively). In the lower 
panel we show the difference between the Sersic-law and core-Sersic 
models. 

is detected, while neither the ELEGs nor the star form- 
ing galaxy 3C 198 have an NIR nucleus. The different 
detection rates in the FR Is and FR lis are not a con- 
sequence of the redshift distribution of the 3CR catalog. 
All FR II sources (HIG, LIG, and BLO) in our sample 
share a similar median redshift, and only the FR Is have 
a lower median z, because of their lower total power. 

Due to the presence of upper li mits, we use the Kaplan - 
Meier product-limit estimator (jKaplan fc Meierl Il958h . 
which provides the mean value for a distribution with 
censored data, where the censoring is random, to mea- 
sure the average value of NIR nuclear luminosities of each 
class (reported in Table [3]). The BLOs have, on average, 
the brightest near infrared nuclear sources in the sam- 
ple in luminosity (<Lir> — 1.6xl0 29 erg s _1 Hz -1 ). 
The FR I nuclei are on average the faintest (<Fir> = 
4.0xl0 26 erg s" 1 Hz" 1 ). FR II LIGs and HIGs share 
a similar range of NIR luminosity and appear to bridge 
the gap between BLOs and FR Is. Undetected NIR nu- 
clei are present for Fir < 10~ 27 erg cm" 2 s" 1 Hz" 1 and 
Fir < 2.5 xlO 28 erg s _1 Hz" 1 . However about one half 
of the detected nuclei have fluxes and luminosities below 
these values. 

4.2. Nuclei in the optical- infrared planes 

It is instructive to compare the NIR properties with 
those of the optical nuclei. In Fig. HI left panel, we plot 
F -F/fl and F -Fir for all objects in the sample. Not 
surprisingly, because of the relative proximity of the two 
bands, most of the points occupy a narrow strip paral- 
lel to the bisectrix of the planes. However, we have to 
be cautious in taking this plot at face value as it is sen- 
sitive to the presence of nuclear absorption and also to 
nuclear variability, since infrared and optical data were 
not obtained simultaneously. 

The FR I galaxies (empty yellow circles) cover the 
whole range in optical and NIR fluxes (~5 dex). The 
FR II LIGs (filled blue circles) are clustered in small op- 
tical and NIR flux range (~2 dex) in the lower-left slide 
of the plane. The HIGs (filled green squares ) are also 
clustered, but with a slight tail towards higher flux, cov- 
ering a wider range in both bands (~3 dex) with respect 
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TABLE 2 

Basic data for all 3CR source considered 
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3C 310 


0.054 


FR I 


B09 


-21.26 


-24.10 


G88 


-28.24 


<-27.07 


<27.69 


39.95 


B09 


cS 


3C 314.1 


0.1197 


ELEG 


B09 


-21.98 


<-26.00 


G88 


<-28.81 


<-27.37 


<28.09 


39.49 


B09 


cS 


3C 315 


0.1083 


FR I 


B09 


-21.75 


-23.75 


G88 




-27.76 


27.52 


40.69 


B09 


cS 


3C 317 


0.0342 


FR I 


B09 


-21.33 


-23.41 


M93 


-27.67 


-27.05 


27.31 


40.25 


B09 


cS 


3C 319 


0.192 


LIG 


B09 


-21.82 


<-25.80 


Z95 


<-28.33 


<-27.77 


<27.96 


<39.91 


B09 


cS 


3C 321 


0.096 


HIG 


B09 


-21.87 


-24.41 


Z95 




-27.64 


27.55 


40.73 


B09 


cS 


3C 323.1 


0.264 


BLO 


B09 


-22.00 


-24.28 


Z95 


-25.83 


-25.77 


30.20 


42.45 


B09 


S 


3C 326 


0.089 


LIG 


B09 


-22.11 


-24.78 


Z95 


<-28.04 


-27.87 


27.25 


40.23 


B09 


cS 


3C 332 


0.1515 


BLO 


B09 


-22.02 


-24.90 


Z95 


-26.76 


-26.08 


29.47 


41.57 


B09 


s 


3C 338 


0.03035 


FR I 


B09 


-21.33 


-23.99 


G88 


-27.78 


-28.17 


26.08 


39.44 


B09 


cS 



Baldi R. D. et al. 



TABLE 2 

Continued 



Name 


z 


FR - E.L. 


Ref 


log FiYSMHz 


log F r 


Ref 


log F 


log Fir 


log Li R 


log L[ 


Ref 


Profile 


3C 346 


0.162 


FR I 


S91 


-21.96 


-23.66 


G88 


-27.42 


-27.01 


28.73 


41.59 


G94 


S 


3C 348 


0.154 


ELEG 


B09 


-20.45 


-25.00 


M93 


-28.88 


<-27.34 


<28.34 


40.17 


B09 


cS 


3C 349 


0.205 


HIG 


J97 


-21.88 


-24.60 


H98b 


-27.62 


-27.39 


28.39 


41.82 


R89 


S 


3C 353 


0.03043 


LIG 


B09 


-20.63 


-23.69 


Z95 


<-28.70 


<-27.99 


<26.24 


40.05 


B09 


cS 


3C 357 


0.167 


LIG 


B09 


-22.09 


-25.15 


Z95 


<-28.73 


<-27.94 


<27.68 


40.70 


B09 


cS 


3C 371 


0.051 


BLLac-LIG 


B09 


-22.43 


-23.04 


P81 


-25.56 


-24.91 


29.75 


40.84 


B09 


S 


3C 379.1 


0.256 


HIG 


B09 


-22.13 


-25.26 


S85 


<-28.74 


<-28.05 


<27.89 


41.51 


B09 


cS 


3C 381 


0.1605 


HIG 


B09 


-21.78 


-25.11 


Z95 


<-28.08 


-27.24 


28.36 


42.12 


B09 


S 


3C 382 


0.0578 


BLO 


B09 


-21.70 


-23.64 


Z95 


-25.05 


-24.91 


29.86 


41.65 


B09 


S 


3C 388 


0.091 


LIG 


B09 


-21.61 


-24.10 


Z95 


-27.88 


-27.83 


27.31 


40.54 


B09 


S 


3C 390.3 


0.0561 


BLO 


B09 


-21.32 


-23.36 


Z95 


-25.70 


-25.42 


29.32 


41.96 


B09 


S 


3C 401 


0.201 


LIG 


B09 


-21.68 


-24.27 


Z95 


-28.15 


-28.26 


27.51 


40.76 


B09 


S 


3C 402 


0.0239 


LIG 


B09 


-22.00 


-24.29 


Z95 


-27.43 


<-27.87 


<26.15 


<39.34 


B09 


cS 


3C 403 


0.059 


HIG 


B09 


-21.75 


-24.92 


Z95 


-28.14 


-26.97 


27.81 


41.62 


B09 


S 


3C 403.1 


0.05540 


LIG 


J97 


-21.87 








<-27.97 


<26.76 






cS 


3C 405 


0.05607 


HIG 


075 


-19.06 


-23.46 


Z95 




-27.84 


26.90 


41.96 


075 


cS 


3C 410 


0.248500 


FR II 


D96 


-21.46 


-22.83 


C08 




-26.25 


29.68 






S 


3C 424 


0.126988 


FR I 


B09 


-21.85 


-24.67 


B92 


<-28.61 


-27.77 


27.64 


40.60 


B09 


S 


3C 430 


0.0556 


LIG 


B09 


-21.47 


-24.77 


S84 




<-27.60 


<27.14 


40.23 


B09 


cS 


3C 433 


0.1016 


HIG 


B09 


-21.25 


-25.24 


G88 




-26.51 


28.73 


41.50 


B09 


S 


3C 436 


0.2145 


HIG 


B09 


-21.75 


-24.61 


Z95 




<-28.06 


<27.76 


41.25 


B09 


cS 


3C 438 


0.29 


LIG 


B09 


-21.35 


-24.60 


Z95 


<-29.18 


<-28.28 


<27.76 


<41.08 


B09 


cS 


3C 442 


0.0263 


FR I 


B09 


-21.79 


-25.68 


G88 


-28.82 


-27.75 


26.35 


39.14 


B09 


cS 


3C 445 


0.057 


BLO 


B09 


-21.61 


-23.42 


Z95 


-25.49 


-25.31 


29.45 


42.39 


B09 


S 


3C 449 


0.01708 


FR I 


B09 


-21.94 


-24.44 


G88 


-27.53 


-27.22 


26.58 


39.13 


B09 


S 


3C 452 


0.0811 


HIG 


B09 


-21.26 


-23.81 


Z95 


-28.16 


<-27.89 


<27.16 


41.19 


B09 


cS 


3C 459 


0.2199 


BLO 


B09 


-21.59 


-22.82 


Z95 




-26.65 


29.18 


41.73 


B09 


S 


3C 465 


0.0301 


FR I 


B09 


-21.42 


-23.57 


G88 


-27.50 


-27.39 


26.86 


39.71 


B09 


cS 


NGC 6251 


0.02471 


FR I 


P84 


-22.00 


-23.05 


J86 


-26.82 


-26.28 


27.80 


39.70 


S81 


cS 



Description of the table: Col.(l): name; Col. (2): redshift; Col. (3): Fanaroff & Riley and optical spectral classifications. For simplicity, 
we mark as FR I all radio galaxies with FR I-like radio morphology and LIG-like optical spectrum, as LIG all radio galaxies with 
FR II-like radio morphology and LIG-like optical spectrum, as HIG (BLO) all radio galaxies with FR H -like radio morphology and 
HIG-Iike optical spe c trum with t he absence ( prese nce) of bro ad emission lines. Re ferenc es (Col. 4): B 09 Buttiglione^t^ay l|2010| ). J97 
IJackson fc Rawlined jl997Tl, H98 a lHardcastTel fl9§8l), Z95 from lZirbel fc Bauml fl995h . S91 ISpencer et al.l 1 119911) , 075 IOsterbrock fe Miller! 
fl975h , D96 lde Koff et all {1996ft , P84 IPerlev et al.l 119841) ; Col.(5): flux density at 178 MHz in erg s" 1 cm -2 Hz -1 (in these logarithmic 
units -23 corresponds to 1 Jy); Col. (6): nuclear radio fl ux density at 5 GHz in erg s — 1 cm" ~ 2 Hz" 1 . Referen ces (C ol. 7): B92IBlack et alj 
J1993) (8.3 GHz) , G88: IGiovannini et al.1 iTM&i. L9 1: ILeahv fc Perlevl Jl99lD . M 93: IMorganti et al l Jl993h . R75 IRilev fc Poolevl 11975ft"T 
T96 ITavlor et all ltl996h. Z95 IZirbel fc Bauml [199511. N98 INilssonl 119981), H9 8b Hardc astliTet al.l 119981) (8.4 GHz) F97: I Fernini et all 
119971). V82: Ivan Breugel & Jager? Tl982 ) . H99: Hardcastlc (1999), B89: Baum & Hcckma3Tl98p), F78: Fomalo mTfc Bridle! 119781), A95 : 
lAkuior fc Garringtonl 119951 ), P81: IPearson fc Readheadl 119811) . S85: ISpangler fc Sakurail 119851) (1.5GHz). S84: ISpangler et al.l jfl984l) . 
J86: IJones et alj 119861 ). C08: A. Capetti, priv. com.. Radio core data at frequencies different from 5 GHz were converted to 5 GHz using 
a flat (a = 0) spectral index; Col. (8): nuclear optical flux density in erg s — 1 cm " 2 Hz -1 IChiaberge et al.l ll999, Chiabcrgc ct al. 2000a] and 
IChiaberge et al.ll2002al) in R band; Col.(9): nuclear 1.6 /im near infrared flux density in erg s 1 cm -' 2 Hz" 1 ; Col. (lO): nuclear 1.6 fim near 
infrared lu minosity in erg s" 1 H z" 1 ; Col.(ll): [O III1AA5007 lum inosi ty in erg s" 1 . Referen ces ( Col.12): B09 IButtiglione et alj 120091 ), 
Q75 IOsterbrock fc Millerl 1T975D , S81 IShuder fc Osterbrockl JT98J), L96 ILawrence et alj 119961 ). R89 IRawlings et alj 119891 ). Col.(13): type 
of model used to fit the NIR surface brightness profile. 
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TABLE 3 
Detections of nuclei in the sample 



Class 


Nuclei 


Upper lim. 


Tot. 


% 


Log hjn 


FR I 


21 


5 


26 


81 


26.7 ± 0.2 


FR II 


42 


31 


73 


58 




LIG 


11 


14 


25 


44 


26.9 ± 0.2 


HIG 


15 


13 


28 


54 


27.5 ± 0.2 


BLO 


15 





15 


100 


29.2 ± 0.1 


ELEG 





3 


3 







SFG 





1 


1 







UNDEF 


1 





1 


100 




BLLac 


1 





1 


100 




Tot. 


64 


36 


100 


64 





Note. — Detection rate of 1.6 £im near infrared nuclei for each 
class. The ELEG (Extremely Low Ionization Galaxies) are 3C 28, 
3C 314.1, and 3C 348. The BL Lac is 3C 371, the Star forming 
Galaxy (SFG) is 3C 198 and the Undefined object (UNDEF) is 
3C 410 because of the lack of spectral data. LIGs, HIGs, BLOs, 
ELEGs, SFG and UNDEF belong to FR II class. Column de- 
scription: (1) class; (2) number of detected nuclei; (3) number of 
undetected nuclei; (4) total number of objects for each class; (5) 
percentage of detection; (6) the mean value and their errors of the 
NIR luminosities for each class with the Kaplan-Meier statistics. 

to FR II LIGs. The BLOs (filled red triangles) show the 
highest optical and NIR fluxes of the sample. 

In the Lo-L/fl plane (Fig. 21 middle panel) we find 
that FRIs exhibit the lowest luminosities. FR II LIGs 
and HIGs are present at increasing optical and NIR lumi- 
nosities. Finally, the BLOs are by far the most powerful 
objects. 

To highlight the presence of outliers in the Lo-^ir 
plane, in Fig. |4] (right panel) we show the ratio between 
NIR and optical nuclear luminosities plotted vs. the NIR 
nuclear luminosity. Most objects lie within ^0.5 dex of 
the dashed line representing i/Lir/i/L = 1. Note that 
most of the nuclei that are found above this line (i.e. 
showing an NIR excess with respect to the optical emis- 
sion) belong to the HIG class. We discuss the origin of 
this behavior of the HIGs in more detail in Sect. 15.41 

4.3. Nuclei in the radio-infrared planes 

In order to investigate the origin of the NIR nuclei, we 
now explore possible correlations between the properties 
of the infrared nuclear emission and those of the radio 
cores. Fig. 03 left panel, shows the IR nuclear flux Fir 
plotted vs. the radio core flux F r . The FR Is (empty yel- 
low circles) cover the whole range of flux spanned by the 
sample, both in the radio and in the NIR. HIGs (filled 
green squares) and FR II LIGs (filled blue circles) are 
clustered in the central region of the plane. The BLOs 
(filled red triangles) are located in the top-right quad- 
rant, at higher radio and NIR fluxes. A similar distri- 
bution of the points is evident in the plane defined by 
the radio and the NIR core luminosities (Fig. [5l right 
panel). However in the L r -L/ij plane the distributions 
are stretched because of the common dependence of the 
two quantities on distance. The FR I galaxies have radio 
core luminosities that are on average lower than those 
of FR lis. This is expected because the Fanaroff- Riley 
morphological radio classification in FR I and FR II cor- 
responds to different total luminosities at 178 MHz, and 
a positive trend links total and core radio luminosity 
(|Giovannini et al.lll988[ ). 

A clear trend is visible for FR Is both in the F r -FiR 



and L r -L]R planes (Fig. [5]). FR II LIGs nuclei also show a 
trend in both planes, even if they cover a smaller range in 
radio core flux (2 dex for the NIR detected nuclei) with 
respect to FR Is. Conversely, there is no evident link 
between radio and NIR emission for HIGs and BLOs. 

In order to assess the presence of correlations, we 
performed a statistical analysis using t he Astronomy 
Survi val Analysis (ASURV) package ()Lavallev et al.l 
119921) . This package is available under IRAF/STSDAS 
and provides us a tool to deal with ce nsored data 
(|Feigelson fc Nelson! Il985l llsobe et al.lll986ft . Note that 
the method we use is based on the assumption that the 
censored data are uniformly distributed. In order to de- 
termine whether this approach is viable, we apply the 
censoring analysis only when the above assumption is 
satisfied. We initially check the location of the non- 
detected data points with respect to the entire range of 
available data. 

We used the "schmittbin" task (|Schmittlll985h to cal- 
culate the associated linear regression coefficients for two 
set of variables. Operatively, we carried out this proce- 
dure twice, obtaining two linear regressions: firstly con- 
sidering the former quantity as the independent variable 
and the latter as the dependent one, secondly switching 
the roles of the variables. The best fit is represented by 
the bisector of t hese two reg r ession lines. This followed 
the suggestion of llsobe et al.l (|1990f ) that consider such a 
method preferable for problems that require a symmet- 
rical treatment of the two variables. In order to esti- 
mate the quality of the linear regression, for small data 
set s (N < 30) we also derived the generalized Kendall's 
r (|Kendalll 1 19831 ) between the two variables, using the 
"bhkmethod" task. Otherwise, we used Spearman's rank 
order correlati on coefficient, using the "spearman" task 
(|Akritaslll9"89l) only for large data sets (N > 30, for the 
cases of HIGs and the whole sample). The statistical 
parameters of all the linear regressions are reported in 
Table H 

The trend shown by FR Is (Fig. [5j empty yellow circle 
points) in the F r -F ir plane corresponds to a linear corre- 
lation coefficient 8 of r = 0.79. The probability to obtain 
this level of correlation from two independent variables 
is P = 2.2 xl0~ 5 . The statistical parameters are even 
better in the radio-infrared luminosity plane, where the 
linear correlation coefficient is r = 0.85, with an asso- 
ciated probability P = 8.9xl0~ 7 . The dispersions of 
these correlations are ~0.51 and ~0.48 dex, for fluxes 
and luminosities, respectively. The generalized Kendall's 
r coefficient, which allows us to include censored data, 
is higher for the luminosities (r = 1.18) than that for 
the fluxes (r = 1.00). The associated probabilities of no 
correlation, are 0.0002 and <0.0001, respectively. Note 
that the slopes of both linear regressions are consistent 
with unity. 

We now turn our attention to the possible presence of 
correlations for FR II galaxies. FR II LIGs show a linear 
correlation in both planes of Fig. [5] Considering only 
the detected nuclei, the linear correlation coefficient is 

8 This is derived using the VLA measurements for the radio 
cores. If we consider the radio cores fluxes from VLBI, we ob- 
tain a similar correlation, statistically indistinguishable from that 
obtained using the VLA cores. This is because the core fluxes 
measured with VLA and VLBI differ only marginally (see the Ap- 
pendix). 
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Fig. 4. — In left panel, logFo-logi 7 /^ (in units of erg s 1 cm 2 Hz 1 ); in median panel, logLo-logLrjj; in right panel, logLm/ L -\ogL[R 
(in units of erg s — 1 Hz ). In these plots the empty yellow circles represent the FR Is, the filled blue circles the FR II LIGs, the filled green 
squares the HIGs and the filled red triangles the BLOs, the purple asterisks the ELEGs, the empty black triangle the BL Lac (3C 371), 
and the filled pink reversed triangle the SFG (3C 198). The solid or dashed lines present in these plots are the line representing the unity 
of the ratio i/hjn/uh . In the right-panel plot, although the same factor hjfj appears in both axes, its aim is not to find a relation or a 
trend, but to better point out the outliers in \og(Lm/ L ). 



r = 0.89 in the ¥ r -Fm plane, with a probability of a 
no correlation of P = 1.0xl0~ 4 . The censoring statis- 
tics gives a generalized Kendall's coefficient r = 0.77 
with a probability, P = 0.0012, that a fortuitous cor- 
relation appears at a level measured by our test statistic. 
In L r -L/fl plane for FR II LIGs a linear relation (dot- 
ted line in Fig. [SJ right panel) is also statistically sup- 
ported: the linear correlation coefficient is r = 0.93 (P = 
1.1 xl0~ 5 ). The radio-infrared core luminosity relation 
for FR II LIGs has a slope of 0.9±0.3. 

Therefore, an important result of our analysis is that 
the correlations found for FR Is and FR II LIGs nuclei 
are consistent with each other within the errors. This 
result is shown in Fig. [6] where for more clarity we have 
marked these two classes with different symbols. 

For HIGs and BLOs the absence of a trend is confirmed 
by the results of the statistical analysis and in particular 
by the high probabilities that their distribution in the 
L,.-L/fl plane is random (see the statistical parameters 
in Tab. SJ| . 

4.4. Spectral indices 

In order to derive the multi-wavelength properties of 
the nuclei, it is useful to analyze the broad-band radio- 
infrared spectral index a r -m. In Fig. [7] we show the em- 
pirical distribution in a r —m for the four different classes 
(FR I, FR II LIG, HIG and BLO). Note that the average 
value of the radio-IR spectral index distributions shifts 
towards lower values from FR II LIGs to the BLOs, with 
FRIs and HIGs having intermediate values. This indi- 
cates that with respect to the IR flux, the FR II LIGs 
are the most radio-powerful objects. Another intriguing 
aspect worth noting is that even though FR II LIGs and 
FR Is behave similarly in the radio-IR plane, as they 
define the same correlations, their distribution of a r -m 
are slightly different. In fact, FR II LIGs have on aver- 
age a slightly higher a r -iR than the FR I galaxies. On 
the other hand, HIGs and BLOs show an near infrared 
excess with respect to the radio emission and therefore 
their a r -m spectral indices are on average lower than 
for the other two classes. 

To assess the presence of any significant differences 



between the spectral index distribution of the different 
classes, we can make use of the 'twosampt' task, avail- 
able in the ASURV package. The task computes several 
nonparametric two-sample tests, giving a variety of ways 
to test whether two censored samples are drawn from the 
same parent population. 

We utilize the Peto-Prentice statistic (|Peto fc Petol 
Il972h . which is a generalization of s tandard test for un- 
censored data, such as t he Wilcoxon (IWilcoxqnl | 19451) and 
Kolmogorov-Smirnov (|Chakravarti fc Rovl I1967D test. 
This statistical test quantifies a distance between the em- 
pirical distribution functions of two samples. The null 
distribution of this statistic is calculated under the null 
hypothesis that the samples are drawn from the same 
distribution. 

We statistically measure the different distributions of 
spectral indices between FR Is and FR II LIGs and be- 
tween HIGs and BLOs. We can assert that FR Is and 
FR II LIGs are not drawn from the same parent popula- 
tion at the 97% confidence level. The same holds for the 
HIG and BLO classes at the 99% confidence level. 

4.5. Comparison with FR Is' radio-optical correlation 

IChiaberge et al.l (|1999f ) measured the optical nuclear 
fluxes for most FR I radio galaxies in the 3CR catalog. 
The optical detection rate is similar to that found for 
the FR Is of our sample. For the sub-sample of objects 
for which the total radio luminosity at 178 MHz is lower 
than the FR I/FR II break (corresponding to \^n%MHz 
~ 2xl0 33 erg s _1 Hz -1 ), they found a tight linear corre- 
lation between the optical and radio nuclear fluxes with 
a dispersion of ^0.4 dex. In order to compare the radio- 
infrared and radio-optical correlations we selected the 15 
FR Is in common with the sample of [Chiabcrg e et al.l 
(1999) (namely 3C 29, 3C 31, 3C 66B, 3C 83.1, 3C 84, 
3C 264, 3C 270, 3C 272.1, 3C 274, 3C 296, 3C 317, 
3C 338, 3C 442, 3C 449, and 3C 465) 9 . We estimate the 
linear regression between radio and near infrared fluxes 

9 Although 3C 277.3 satisfies the selection criteria, we excluded 
the object from this sample since from an careful inspection of its 
radio maps we revised its radio-morphological classification from 
FR I to FR II. 
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Fig. 5. — Left panel: radio core flux vs. near infrared nuclear flux (in units of erg s 1 cm 2 Hz 1 ). Right panel: radio core luminosity 
vs. near infrared nuclear luminosity (in units of erg s — 1 Hz -1 ). In these plots the empty yellow circles represent the FR Is, the filled blue 
circles the FR II LIGs, the filled green squares the HIGs and the filled red triangles the BLOs, the purple asterisks the ELEGs, the empty 
black triangle the BL Lac (3C 371), the filled pink reversed triangle the SFG (3C 198), and the filled red diamond is the undefined object 
(3C 410). The solid line represents the linear correlation found for FR Is and FR II LIGs. 
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Fig. 6. — Radio core luminosity vs. near infrared nuclear lumi- 
nosity (in units of erg s _1 Hz -1 ) for FR Is (empty yellow points) 
and FR II LIGs (blue points). We show the linear correlations: the 
solid line is obtained for FR Is and FR II LIGs, the dotted line for 
FR Is only, and the dashed line for FR II LIGs only. 

for this sub-sample: the correlation coefficient is r = 0.82 
(P = 2.2 xl(T 4 ) and the dispersion is ~0.49 dex. We 
tested that the dispersions of the two relations, F r -Fm 
and F r -F„, are not statistically different at a significance 
level of 95%. 

5. DISCUSSION 

In the following we consider the properties of each 
group of 3CR radio galaxies separately. 



5.1. FR I galaxies 

The presence of a linear correlation between FR I ra- 
dio and NIR cores, both in flux and luminosity, argues 
for a real physical connection between the two variables. 
This is indeed not the result of a common dependence 
on redshift of the two quantities, as shown by the fact 
that the statistical parameters of b o th cor relations are 
similar. Following IChiaberge et al.l ()1999l ). that inter- 
preted the optical-radio core correlation as the indica- 
tion that synchrotron emission dominates the nuclear 
emission at both wavelengths, we argue that the corre- 
lations between the radio and the NIR nuclei described 
here can be explained by the same physical mechanism. 
More precisely, the presence of such a tight correlation 
between the NIR and radio cores indicates that non ther- 
mal (synchrotron) radiation dominates their emission at 
all considered wavelengths, i ncluding the X-ray as shown 
bv lBalmaverde et al.l (|2006[ ). iCapetti et all (|2007| ) found 
that the optical HST nuclei of 9 3CR/FR I, which also be- 
longs to our sample, are highly polarized. This strength- 
ens the synchrotron origin of their nuclear emission. 

The correlations among NIR, VLA and VLBI cores 
and the fact that the fluxes of the VLA and VLBI cores 
(see Appendix) are essentially equal suggest that the NIR 
cores are typically produced on scales consistent with or 
even smaller than the size of the VLBI cores. In fact, ac - 
cording to various jet models fe.g.. lGhisellini et"aill985h . 
it is likely that the radio and the IR/Optical radiation 
is produced in the innermost region of the jet, close to 
the black hole (< 10 16 cm), while the radio is emitted on 
larger scales. 

The dominance of synchrotron radiation even in the 
NIR indirectly confirms that accretion in most FR I nu- 
clei occurs at low rates and/or with low radiative ef- 
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TABLE 4 

Summary of the statistical analysis of the linear regressions. 



Relation 


Class 


N d 


N c 


rd 


Pr„ 


Pd or r d 


Pp. or P TH 


p or t 


P P or P T 


rms x 


rmsj, 




m 


q 


Fr-F/fl 


FR I 


21 


5 


0.788 


2.2xl0" 5 


1.181 


0.0002 


0.997 


0.0002 


0.477 


0.543 


1 


0±0.2 


-3±5 


Fr-F/fl; 


FR I a 


15 





0.834 


0.00011 


1.238 


0.0013 


/ 


/ 


0.442 


0.542 





8±0.2 


-8±6 


Fr-F 7fl 


FR I b 


15 


1 


0.855 


4.9xl0" 5 


1.257 


0.0011 


1.233 


0.0006 


0.394 


0.407 


1 


0±0.3 


-3±6 


Fvlhi-FlR 


FR I 


15 


1 


0.866 


3.0xl0~ 5 


1.239 


0.0012 


1.217 


0.0007 


0.380 


0.291 


1 


3±0.3 


5±5 


Fr-Fjfl 


LIG 


11 


12 


0.890 


0.00010 


1.490 


0.0014 


0.775 


0.0012 


0.168 


0.298 





6±0.3 


-14±15 


Fr-Fj fl 


FR I+LIG 


32 


17 


0.781 


4.1xl0~ 8 


0.743 


<0.0001 


0.645 


<0.0001 


0.445 


0.511 





9±0.2 


-7±14 


Fr-F Jfl 


HIG 


15 


12 


0.150 


0.59 


-0.095 


0.804 


0.296 


0.218 


0.542 


0.544 


1 


2±0.4 


3±24 


Fr-Fj fl 


BLO 


15 





0.348 


0.20 


0.514 


0.181 


/ 


/ 


0.680 


0.727 


1 


0±0.4 


-3±9 




FR I 


21 


5 


0.854 


8.9xl0" 7 


1.324 


<0.0001 


1.077 


<0.0001 


0.420 


0.532 


1 


0±0.2 


-3±6 




FR l b 


15 


1 


0.898 


5.4xl0" 6 


1.276 


0.0009 


1.167 


0.0008 


0.392 


0.417 


1 


0±0.2 


-3±5 




FR I 


15 


1 


0.903 


4.1xl0" 6 


1.333 


0.0005 


1.183 


0.0007 


0.384 


0.347 


1 


1±0.2 


-7±5 


Lr-L/ii 


LIG 


11 


12 


0.931 


l.lxlO -5 


1.491 


0.0014 


0.696 


0.0011 


0.211 


0.306 





9±0.3 


0±12 


Lr-L/ii 


FR I+LIG 


32 


17 


0.890 


1.8xl0~ 8 


0.900 


<0.0001 


0.688 


<0.0001 


0.363 


0.487 





8±0.2 


1±18 


Lr-L/ii 


HIG 


15 


12 


0.503 


0.056 


0.514 


0.181 


0.461 


0.054 


0.637 


0.606 


1 


3±0.5 


-13±27 


Li-L/ii 


BLO 


15 





0.2398 


0.39 


0.152 


0.692 


/ 


/ 


0.488 


0.802 





8±0.5 


4±14 



Note. — Column description: (1) relation studied; (2) class of radio-galaxies considered for the relation; (3) number of detected 
nuclei included in the relation considering both the quantities; (4) number of censored nuclei included in the relation considering both the 
quantities; (5) r d linear correlation coefficient for the relation obtained only with the detected data; (6) P r<J probability associated with r d 
that the correlation is not present; (7) Pd Spearman's rank coefficient (if > 30) or r^, generalized Kendall's tau (if N,j < 30), for the 
relation obtained only with the detected data; (8) P Pd (or P Td ) probability associated with p d (or r^) that the correlation is not present; 
(9) p Spearman's rank coefficient (if N^+c > 30) or r, generalized Kendall's tau (if Nj +C < 30), considering all the data; (10) P p (or P T ) 
probability associated with p (or r) that the correlation is not present; (11) rmsj, rms on x-axis for the relation obtained only with the 
detected data; (12) rms y , rms on y-axis for the relation obtained only with the detected data; (13)-(14) slope coefficient m and intercept 
coefficient q for the linear regression (y = m*x + q) obtained considering all the data. In the rows marked with a the considered sample 
is composed of the FR I objects with \^\i%mHz < 2xl0 26 W Hz -1 shared with sample studied by Chiabergc ct ah| 1119991) . In the rows 
marked with b the considered sample is composed of the FR Is whose the VLBI radio core measurement is available. With LIG we consider 
FR II LIG. 
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ficiency (e.g.. iBaum et al.lll995l IChiaberge et al] |1999l 



I Allen et alJ I2006L iBalmaverde et all 120081 ) . However, 
there are exceptions to this result: the clearest exam- 
ple is possibly 3C 120 (a radio galaxy which is for- 
mally not part of 3CR sample) which d isplays a ra- 
dio m orphology typical of FR Is (e.g. IWalker et alJ 
|2001[ ) while its nuclear emission shows the clear signa- 
ture of broad lines and rad iatively efficient accretion 
(|Phillips fe Osterbrockl 119751) . On the other hand, ex- 
ceptions are also present among BL Lacs, the putative 
"beamed" counterpart of LIG population. A number of 
BL Lacs also s how (faint) narrow an d/or broad emis- 
sion lines (e.g. IVermeulen et al.lll995l ) and others have 
FR Il-like extende d radio emission (e.g. lAntonuccllll986l : 
iKharb et al.ll2010h 

IChiaberge et alJ ()1999l ) pointed out that the dispersion 
of the radio-optical correlation (~ 0.4 dex) is consistent 
with dust absorption, provided that the Ay is randomly 
distributed among the sources and does not exceed ~ 2 
mag. Since the extinction in H-band is 4.5 times lower 
than in the optical R-band, the observations presented 
in this work are far less sensitive to the presence of dust 
along the line of sight than those in IChiaberge et alJ 
(1999). Therefore, one would expect the dispersion of 
the infrared-radio correlation to be smaller than that of 
the radio-optical, if that is due to absorption. Our re- 
sults instead show that the rms of the two correlations 
are similar. Thus, the most likely explanation for the 
scatter is that variability of the nuclei plays a fundamen- 
tal role, while the effect of absorption is less significant. 
Indeed, this is supported b y the variability observe d for 
the nuclear flux of 3C 317 (IChiaberge et alj|2002bl) an d 
of 3C 274 (jTsvetanov et al l 119981 iPerlman et all 120031 ) . 
Conversely, the presence of dust appea rs to play a role 
in dim ming the nuclear flux on the UV fChiabcrg e et all 
2002b). Note that data in the literature indicate that 
the radio core variabili ty is typically at 10% over a time 
scale of several years (iHine fe Scheuedll980l: lEkers et all 
119831 ). Therefore the main source of the scatter in the 
radio-NIR correlation is mostly likely due to a variable 
NIR nuclear component. 

Furthermore, the region of the jet base that produces 
the optical/NIR emission in FR Is is likely to be mildl y 
relativistic (e.g.. IChiaberge et allfl999t IXu et al.lll999bD . 
Therefore, we consider the effects of Doppler boosting on 
our results. Firstly, the range in NIR nuclear luminosity 
might be enhanced by relativistic beaming, which would 
boost (or de-boost) the emission of FR I nuclei, depend- 
ing on the viewing angle. Objects seen under a small 
angle between the line-of-sight and the jet axis would ap- 
pear boosted, while objects observed closer to the plane 
of the sky would be de-boosted. Secondly, nuclear vari- 
ability might be enhanced by relativistic beaming when 
the flux variation is associated with 1) a change in the jet 
bulk Lorentz factor, or 2) a change in the jet direction, 
which would result in a different angle between the jet 
direction and the line-o f-sight. 

However, as noted bv lChiaberge et aD(|1999l ). the bulk 
Lorentz factor of the component responsible for the nu- 
clear optical/NIR emission is most likely around r ~ 2 
and the source is possibly a "slower" shear layer in the 
jet. Therefore, the effects of relativistic beaming should 
not dramatically affect our results. 

Are the NIR data telling us something new about the 



presence of optically and geometrically thick obscuring 
tori in FR Is with respect to the optical observations? 
The results derived considering the UV, optical, and NIR 
bands are only apparently in contradiction. In fact, the 
median absorption toward FR I nuclei derived from UV 
observations is Ay ~ 1.3. This corresponds to only 0.2 
mag in H band, negligible with respect to a dispersion of 
~0.5 dex, and to ~ 3 mag at 2500 A. 

The detection rate of FR I nuclei in t he NIR (81%) is 
simila r to that found in the R band by IChiaberge et all 
(fl999l) . All sources with an optical nucleus also show a 
NIR counterpart, except 3C 310. Instead, in 3C 75N a 
faint dust lane prevents the study of its nuclear regions in 
the optical, while its NIR core is clearly detected. Fur- 
thermore, the absence of FR I outliers in the Lm/Lo 
ratio (see right panel of Fig. 2J can be interpreted as a 
hint that the detected FR I nuclei are not affected by 
significant absorption (Ay < 2 mag). 

The reason for the non-detection of some of the FR I 
NIR nuclei might instead be due to a small contrast be- 
tween the nucleus and the host galaxy stellar emission. 
Assuming that the radio and the near infrared emission 
of FR I nuclei are intrinsically tightly correlated, the in- 
tensity of their radio cores provides a robust prediction 
of their infrared nuclear flux. Following the a nalysis on 
the B2 radio-galaxies by ICapetti et all ()2002l ). we plot 
the radio core fluxes of all FR Is against the galaxies 
central surface brightness (Fig. [5J . Four out of the five 
undetected nuclei are indeed among those with the lowest 
predicted contrast against the galaxy emission. 3C 310 is 
the only object that does not have an NIR detected nu- 
cleus and apparently shows an rather larger contrast with 
galaxy emission than that of the other four non-detected 
nuclei. However, since its optical nucleus is detected, its 
NIR non-detection is likely not due to a obscuration. 

Summarizing, we argue that in most cases we have a 
direct view of FR I nuclei and that the few undetected 
NIR central sources are below the detection threshold 
set by the contrast with the underlying host galaxy emis- 
sion. Note that the detection of massive du sty structures 
in the central regions of some FR Is (e.g. lOkuda et all 
l2005llDas et al.ll2005l ) does not contrast with our result 
that the NIR nuclei are mostly unobscured. In fact dust 
is distributed in large kpc-scale geometrically thin disks 
which do not appear to 'hide' the nuclei to our line-of- 
sight. 

5.2. FR II: Low Ionization Galaxies LIG 



IChiaberge et all (|2002aD found that the nuclei of FR II 
LIGs lie on the FRI radio-optical correlation. The sam- 
ple of FR II LIGs available in this work is substantially 
larger than the sample those authors considered, both 
because of the increasing number of FR II LIG identifi- 
cations and of the larger sample size. This allows us to 
establish the presence of a correlation between the radio 
and the NIR emission within the sample of FR II LIGs 
itself. An important result is that the correlations for 
FR Is and FR II LIGs are statistically indistinguishable. 
This strengthens the hypothesis that the LIGs are FR II 
galaxies with FR I-like nuclei. Such a nuclear similar- 
ity implies that non-thermal synchrotron radiation (e.g., 
from a relativistic jet) dominates their nuclear emission 
in NIR band and that accretion may occur with a radia- 
tively inefficient process and/or at low accretion rates in 
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Fig. 8. — Radio core flux density, log F ra[ j; (erg s — 1 cm -2 
Hz -1 ), versus central (on a circular region of r=0''15) surface 
brightness of the host galaxy, log F ga i (erg s — 1 cm -2 A -1 
arcsec -2 ) for all FR I. Objects with detected near infrared nu- 
clei are marked with empty yellow circles, while those with non- 
detection by filled yellow circles. The dotted line marks a constant 
ratio between the radio core and galaxy emission. 

FR II LIGs as well. 

The large fraction of FR II LIGs in our sample (^33% 
of the FR lis) proves that these objects are not rare ex- 
ceptions in the FR II class. This class of objects show 
nuclear properties similar to those of FR Is, whereas their 
radio morphologies span the range from classical double 
to X-shaped, and from com pact source to Fat Double 
(FD). The FD morphology (|Owen fe Laingl 119891 ) con- 
sists of a FR II radio structure with elongated and dif- 
fuse lobes and sometimes with prominent jets, and may 
be considered as an intermediate class between FR Is and 
FR lis. At least half of the FR II LIG sample show a 
FD radio morphology. The large variety of extended ra- 
dio morphologies associated with a central LIG-like AGN 
makes it difficult to understand how the properties of 
the core and extended radio structures are related (e.g. 
!Bicknelll[l98l[T99llGoDal-Krishna fc Wiitall2001h . The 
existence of FR II LIGs as a class shows that radiatively 
inefficient accretion disks are able to produce jets that 
give rise to the typical large scale FR II morphology. 
Another possibility is that FR II LIGs constitute an evo- 
lutionary (or transient) state of the life of the FR II class. 
However, the environment of FR II LIGs is different from 
that of other classes of FR lis, as FR II L IGs are predomi- 
nantly found in dense gr oups or clusters (Chiabcrg e et all 
l2000at lHardcastldl2004l ) . That seems to rule out the evo- 
lutionary scenario. 

The lower detection rate of FR II LIG nuclei with re- 
spect to that of FR Is is most likely due to the reduced 
contrast between their faint IR nuclei and the surround- 
ing stellar emission of the host galaxy, because of the 
higher redshift of FR II LIGs as compared to that of 
FR Is. In fact, the fraction of detected FR II LIG nuclei 



decreases as redshift increases, from 55 % for z<0.15 to 
29 % for z>0.15. However, at this stage we cannot rule 
out that the presence of significant obscuration along the 
line of sight to the FR II LIGs with non-detected nuclei 
might also contribute to the reduce the fraction of de- 
tected FR II LIG nuclei. The best way to probe obscu- 
ration in FR II LIGs is through X-ray spectral analysis. 
FR II LIGs indeed seem to lack s ignificant absorbing co l- 
umn densities, N H ~ 10 23 cm" 2 (jHardcastle et al.l l2006) . 

The fact that the central engines of FR II LIGs and 
FR Is are similar (i.e., they are both unobscured and lack 
a significant broad line region) has important bearings for 
the un ification schemes. As pointed out bv lLaing et al.l 
(1994), the broad distribution in core-dominance is con- 
sistent with FR II LIGs being a randomly oriented sam- 
ple. Since FR II LIGs also appear to lack significant 
broad emission lines, when observed along the jet axis, 
thes e objects would appear as BL Lacs, in agreement 
withlChiabcrg e""et al.l ()2000bl ) in light of their findings in 
the optical band. Such a unification scenario, which was 
previously explored by iJackson fc Walll (|1999l) based on 
the optical spectral properties of FR II LIGs, is in appar- 
ent contrast with zeroth order unification scheme which 
associates all FR II with quasars. However, besides be- 
ing in agreement with our results, it may also account 
for the presence of BL Lacs with an FR II morphology. 

Nevertheless, FR II LIGs do show some differences 
with respect to FR Is. In fact, as shown in Fig. [71 FR II 
LIGs and FR Is have slightly different radio-IR broad- 
band spectral indices. If FR II LIGs are indeed the coun- 
terparts of BL Lacs with FR II radio morphology, the 
spectral index discrepancy can be explained in the frame- 
work of the so-called "blazar sequence" (jFossati et al.l 
HflflSh . Such a scenario predicts a smooth transition 
in the properties of blazars' spectral energy distribution 
(SED). In blazars, the synchrotron emission peak is ob- 
served to shift towards lower frequencies as the radio lu- 
minosity at 5 GHz increases. The luminosity of the syn- 
chrotron emission peak also increases following the same 
trend. This gives rise to a sequence having high energy 
peaked BL Lacs at the lowest luminosities, low energy 
peaked BL Lacs at intermediate luminosities, and flat 
spectrum radio quasars at the highest bolometric powers. 
In Fig. [9] we schematically show what the effect of such 
a scenario on the broad band spectral indices would be. 
As the peak in the SED shifts towards lower frequencies, 
otr-iR steepens. Therefore, the steeper spectral indices 
we found in FR II LIGs might be explained in terms of 
a slightly higher 5GHz core luminosity with respect to 
that of FR Is. FR II LIGs would thus correspond to in- 
termediate power blazars such as the low-energy peaked 
BL Lacs. 

The possibility that, conversely, the change in a r -iR 
is due to the redshift of the SED for the more distant 
objects can be exclud ed using the formula presented by 
iTrussoni et al.l (|2003[ ). The formula, originally derived 
for estimating the effect of relativistic beaming, can be 
used to to account for the effect of redshift on the broad- 
band spectral indices, by substituting R^, the ratio of 
the beaming factors, with (1 + z). The variation in the 
radio-infrared spectral due to the frequency shift for an 
object at z^O moved to z=0.3 is only Aa r ^m ~ 0.02, a 
correction that has a negligible effect on the results. 
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Fig. 9. — The Spectral energy distributions (S EP) of three 
blazars in order to explain the so-called scenario (Fossati ct al. 
1998). In the three blazars, the synchrotron emission peak shifts to- 
wards lower frequencies as the radio luminosity at 5 GHz increases, 
which also corresponds to a higher luminosity of the synchrotron 
emission peak. The dashed lines represent the radio-infrared spec- 
tral indeces of each SED. As the peak in the SED shifts towards 
lower frequencies, a r —jn steepens. 

5.3. FR II: Broad Line Objects BLO 

According to the standard unification scheme, the 
BLOs are objects where the nuclear region i s un- 
obscured along our line o f sight (e.g. iBarthell 119891 : 
lAntonucci fe Ulvestadl lT985). This idea is supported by 
the fact that all the NIR nuclei are detected for BLOs. 

The BLOs show a strong NIR excess, of up to two 
orders of magnitude in luminosity, with respect to the 
correlation defined by the FR I (and FR II LIG) radio- 
infrared nuclei (Fig. [5J right panel). A s imilar excess was 
already found in the optical band by IChiaberge et al.l 
(2002a) interpreted that as thermal emission from the 
accretion disk. 

Let us now explore the origin of their NIR nuclear light. 
The infrared-optical spectral index for BLOs covers the 
range am~o ^0-2 (with the exception of 3C 111 for 
which aiR- = 3.3). Such a large range might result 
fr om a wide distrib ution of spectral indices (see Fig. 11 
in lElvis et al.l H.994). In addition, variability may smear 
the intrinsic distribution. Therefore we consider the me- 
dian value aiR- = 0.95. This value indicates that the 
Spectral Energy Distribution (SED) of BLOs is essen- 
tially flat in a i^L„ representation. 

To better investigate the origin of NIR emission in the 
BLOs of our sample we compare their pro perties with 
those of the radio- loud QSOs (RLQSO) from lElvis et al.1 
(1994). As demonstrated by those authors, the averaged 
SED of QSOs shows a minimum at ^1/im. For wave- 
lengths shorter than 1/xm the light is likely to be dom- 
inated by the accretion disk (i.e., the big blue bump), 
while for longer wavelengths a broad peak is observed, 
which is interpreted as emission from hot dust. The flat 
spectral index of BLOs is an indication that, on aver- 
age, the optical and NIR HST observations are located 
on opposite sides of the SED minimum. Therefore, by 
analogy with RLQSOs, we then conclude that in most 
3CR/BLOs the H band light is mainly produced by the 
high temperature tail of the hot dust emission, located 
at the inner face of the obscuring material. 



We estimate a spectral index of am- — 0.35 from the 
averaged SED of RLQSOs after adopting the mean red- 
shift of BLOs (z = 0.14). The difference of the spectral 
indices between RLQSOs and BLOs can be accounted by 
an NIR/optical flux ratio of the BLOs larger by a fac- 
tor ~2 than that of RLQSOs. This is an indication of 
a relatively larger contribution from the dust in infrared 
emission of BLOs than that observed in RLQSOs. Fur- 
thermore, reminding the reader that the RLQSO have 
optical and NIR luminosities a factor of 100-1000 higher 
than those of the 3CR/BLOs, there is a dependence of 
the optical-infrared spectral index on the luminosity of 
the object. 

5.4. FR II: High Ionization Galaxies BIG 

At odds with the results found for FR II LIGs and 
FR Is, the NIR emission of HIG nuclei does not correlate 
with the radio cores. Furthermore, most of the detected 
nuclei exceeds the emission of FR Is and FR II LIGs at 
a given radio-core luminosity, suggesting an additional 
component dominating over the jet emission, similar to 
what we found for the BLOs. 

However, HIG nuclei are substantially fainter than 
those of the BLOs, on average by a factor of 50. This 
is expected in a standard unification scheme, since the 
HIGs are considered to be the obscured counter parts of 
BLOs. This result indicates that substantial obscuration 
is still present in the NIR, despite the reduced effects of 
dust in the H band. 

In order to investigate the origin of NIR emis- 
sion in HIGs, we fol low the method described by 
iMarchesini et al.l (|2005| ) who studied a complete sub- 
sample of the 3CR radio catalog, which includes all HIGs 
with z < 0.3, observed in K' band from the ground. They 
first estimated the ionizing luminosity, L; orl = \jnlr 
C _1 =0.0675 x ~L[oni] (e.g. total luminosity in nar- 
row lines, Lnlb., estimate d from [O II] and [O III], 
iRawlings fc Saunders! Il99l for HIGs and BLOs assum- 
ing the same covering fa ctor for the two cl a sses, C = 0.01. 
Then, for each source, IMarchesini et al.l (|2005l ) plotted 
Li on /Ljn against L /Ljb, which are both sensitive to 
nuclear obscuration. They found that HIGs do not fol- 
low the 'reddening line' (see caption of Fig. ITUl) in such a 
plane, as one might expect if they were simply "partially 
obscured" BLO nuclei. This lead them to conclude that 
a substantial role i s played by scattered nu clear light. 

With respect to IMarchesini et al.l (|2005l ) our analysis 
can rely on direct measurements of the NIR nuclei and 
these are available for a larger number of HIG, due to 
better sensitivity of the HST images to faint IR nuclei. 
Furthermore, we also have NIR data for the nuclei of 
BLO, thus providing a more robust estimate of the loca- 
tion of unreddened sources. 

In agreement with the results IMarchesini et al.l (2005) 
we found that only a few HIGs are located along the 
reddening line in Fig. [TU] but this is clearly not the case 
for the whole population. 

We then include the effects of scattered nuclear light. 
In this scenario a dusty torus attenuates the direct nu- 
clear light by an extinction Ay :torMS and a fraction / of 
the total nuclear light Ljr is scattered into o ur line of 
sight (see Eq. 3 and in Marchc sini et al.ll2005| ). The ef- 
fects of varying Ay :torMS is as follows: an increase of the 
nuclear obscuration simply corresponds to an increase of 
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Fig. 10. — Ratio between the ionization luminosity L; on and 
z^L/i} versus the ratio between L and L/^. The BLOs are the 
red triangles and the HIGs the green squares. The solid oblique 
line is the 'reddening line' that purely obscured objects should 
follow. The origin of this line is determined by the ratios measured 
for BLOs, i.e. \og(L ion /uh IR ) = 1.17±0.15 and log(L /L/ fl ) = 
-0.28±0.11. The dashed region reproduces the uncertainty of its 
origin. The dotted vertical line corresponds to the average value 
of log(Lo/L/fl) measured for BLOs. The grid of dotted curves 
represent the effects of varying th e nuclear A v ,t and the fr action of 
scattered light / (see the Eq.3 in Marchcsini ct al. (2005)). 

hi on /i>LiiR. The behavior of L /L/a is more complex: 
initially, an increase of the absorption moves the point 
(starting from the average location of BLO) to the right, 
along the reddening line. At larger values of Ay torus the 
direct (transmitted) optical light becomes fainter than 
the scattered component, while the infrared nucleus is 
still substantially unaltered. At even higher values of 
Ay,torus, the infrared transmitted component starts to 
decrease; since the scattered optical component remains 
constant this causes an increase of Lo/L/^, producing a 
turnover in the source path. The track ends at the same 
ratio of L /L/fl found for BLOs, but shifted downward by 
a factor linked to Ay. torus . The location of the turnover 
instead depends on the fraction of scattered light /, mov- 
ing to smaller L /Ljh and to larger L ion /vLiR with de- 
creasing /. 

Because both the IR and the optical emission is 
most likely variable, therefore non-simultaneous mea- 
surements may lead to large errors in the estimate of 
the intrinsic Lo/L/^. Nevertheless, this does not com- 
promise the analysis: while a handful HIGs can be inter- 
preted with being purely obscured BLOs, most of them 
(approximately 3/4 of the sample) require obscuration as 
well as scattered light, in the range of 10<Ay itorus <30 
and 0.02%< / < 10% respectively. This agrees with the 
general results of lMarchesini et all (|2005t ) (they derived a 
fraction of ~70%) and also considering HIGs on a object- 
by-object basis. 

Our analysis of the brightness profiles of the 3CR 
sources always assumed that the nuclear NIR compo- 



nents are unresolved. While for the other classes of 
3CR galaxies the nuclear sources are directly associated 
with the AGN (and thus effectively point-like at our 
resolution), for HIGs our results suggest that they are 
in general dominated by scattered light, a spatially re- 
solved component. This is confirmed by several observa- 
tions of imaging p olarimetry and spectropolarimetry of 
radio -galaxies (e.g. iTadhunter et al.l l2000: Ra mirez et al.l 
2009). We then re-modeled the observed NIR brightness 
profiles taking into account the possibility that the cen- 
tral HIG sources are actually resolved and can be de- 
scribed by a gaussian light distribution. We considered 
in detail the three objects (3C 234 10 , 3C 300, and 3C 349) 
that should have the larger scattered light fraction, be- 
ing located closer to the grid origin in Fig. [lOj For sizes 
larger than ~ O'.'l significant residuals start to appear in 
the model fitting and we consider this value as an upper 
limit to the size of the scattering region. This corre- 
sponds to a range of 300-400 pc for these three objects, 
consistent with the size of the Seyfert polar scattering 
zone. This does not mean that there is no scattered ligh t 
outside this region (as in e.g. 3C 321. lHurt et aH H999V 
but it must be taken as an indication that the bulk of 
scattering occurs within this scale. The presence of low 
surface brightness scattered components (in contrast to 
the host galaxy) can only be probed with imaging po- 
larimetry. 

We conclude that HIG, unlike FR II LIGs, arc con- 
sistent with the interpretation that they host a hidden 
quasar, but they are not simply partially obscured BLO. 
Scattered light must provide a substantial contribution 
to their nuclei in both the optical and NIR bands. 

6. SUMMARY AND CONCLUSIONS 

We have analyzed 1.6 /Ltm near infrared images of 100 
3CR radio galaxies from HST NICMOS in the F160W 
band with z < 0.3, ~90% of the whole 3CR sample. On 
the basis of the radio morphology we divided the sample 
between FR I and FR II and from the point of view of 
their optical spectra we also classified them as Low Ion- 
ization Galaxies (LIG), High Ionization Galaxies (HIG), 
and Broad Line Objects (BLO). 

We measured the nuclear NIR luminosity of all the 
objects, by subtracting the 1-D radial brightness profile 
(modeled with a Sersic or core-Sersic law representing 
the galaxy emission) from the total profile derived from 
the isophote fitting of the images. 

The NIR nuclear emission is found to be strongly cor- 
related with radio cores for FR Is. The linear correlation 
extends over ~ 4 dex in luminosity with a rms of 0.4 dex 
and it has slope consistent with unity, as already found 
for the relations with optical and X-ray nuclei. This cor- 
relation provides further support to the identification of 
their nuclear emission as the synchrotron radiation pro- 
duced in the inner part of the relativistic jet. In turn this 
implies a low contribution from thermal emission, an in- 
dication for the presence of radiatively inefficient disks. 
The similar dispersions for the correlations in NIR and 
in optical with the radio cores ascribes to the nuclear 
variability the dominant role of this scatter, rather than 
to obscuration. 

10 3C 234 is a clear well-known cas e of optica l ly obscured 
AGN with broad Ho in polarized flux HAntonuccil 11981 IT9841 : 
I Young et id]|l998ft . 
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The large number of FR II LIGs present in this sample 
made possible for the first time a detailed analysis their 
nuclear properties in radio-infrared plane. We found a 
strong linear correlation between radio and NIR emis- 
sion. The multiwavelength properties of FR II LIG nu- 
clei are statistically indistinguishable from those of FR Is, 
indicating a common structure of the central engine, de- 
spite the difference in radio morphology and radio power. 
The similar nuclear properties suggest that the nuclei of 
FR II LIGs belong to the same class of those of FR Is, 
i.e., the nucleus is dominated by non-thermal emission 
from a relativistic jet. 

All BLOs have detected NIR nuclei and show an NIR 
excess with respect to the FR I-FR II LIG radio-infrared 
correlation. This excess does not correlate with the radio 
cores, suggesting the presence of an additional infrared 
emitting component. In analogy with the properties of 
radio-quiet QSO, their optical nuclei most likely originate 
from the thermal emission of the accretion disk. Since 
their optical-infrared spectral indices are essentially flat, 
this suggests that thermal infrared emission from hot 
dust dominates in the NIR. 

The origin of the NIR nuclei of HIGs is more complex. 
Although a few HIG nuclei lie on the FR I-FR II LIG 
radio-infrared correlation, most of them show a large NIR 
excess which does not correlate with the radio cores. This 
means that synchrotron emission from the jet base is not 
sufficient to explain their nuclear luminosities. Nonethe- 
less, the HIG NIR nuclei are less luminous than those 
of BLOs, by a factor of 50 on average, suggesting par- 
tial obscuration, as expected by the unified model. Yet 
the infrared-optical luminosity ratios for HIGs are not 
consistent with a model of HIGs as the purely-obscured 
counterpart of the BLOs. We found that this inconsis- 
tency can be solved by including a significant contribu- 
tion from nuclear light reflected in a compact circum- 
nuclear scattering region. The sources which show the 
largest scattered light contribution, have a scatter region 
of 300-400 pc (~ CK'l). 

The analysis of the properties of the NIR nuclear 
sources of 3CR radio-galaxies presented here broadly 
confirms the indications derived from the study of the 
optical nuclei. In particular, the presence of two classes 
of radio-loud AGN, LIG and HIG, with different nuclear 
properties is strengthened. 

FR I radio-galaxies as well as the sub-population of 
sources with FR II morphology characterized by a low 
ionization optical spectrum belong to the same class from 
the nuclear point of view. In these sources the mecha- 
nism of emission related to the presence of non-thermal 
plasma from their radio-jets dominates over the contri- 
bution from thermal radiation related to the accretion 
process. Moving from optical to infrared observations 
has no substantial effect despite the reduced absorption 
at longer observing wavelengths. However the NIR band 
is not the optimal band to investigate the issue of ob- 
scuring tori. The data preseneted here leave the re- 
sults based on optical observations unchanged, bearing 



in mind that for most of FR Is we have a clear view to 
thei r optical/infrared nu cleus. While MIR Spitzer spec- 
tra (jLeipski et al.l [20091 obtained with a large slit size 
(10.7") do show an infrared excess with respect to the 
jet emission in three FR I sources (and in one FR II 
LIG), our results on 3CR sample do not directly require 
the presence of dusty torus on th e line of sight to the F R I 
nucleus. Furthermore, although iLeipski et al.l (|2009f ) ar- 
gue that the dust emission is powered the nucleus, the 
comparison between MIR spectra and HST data is ex- 
tremely difficult because of the different spatial resolu- 
tion of the two instruments. On the other hand, at this 
stage we cannot rule out the presence of a dusty torus 
in these objects because of the higher distance of the 
FR II LIGs than that of FR Is. Despite the lack of such 
information, the nuclear properties of FR II LIG nuclei 
are found to be statistically indistinguishable from those 
of FR Is. How all the sources of the FR I/FR II LIG 
group can be associated with such a large range of radio 
power, while still maintaining their disks in radiatively 
inefficient state, and what causes the transition to the 
FR II morphology in LIGs remain open questions. 

The NIR data provide new insights also for the sec- 
ond class of objects, formed by HIGs and BLOs. The 
BLO NIR nuclei are likely to be dominated by the emis- 
sion from hot dust located in the immediate vicinity of 
the active nucleus. In HIG we found many examples of 
nuclei in which the NIR emission is substantially larger 
than in the optical. Due to the proximity of the two 
bands, the only mechanism that can explain this effect is 
the presence of nuclear structures covered by absorbing 
material, optically thin at least in the infrared. Scat- 
tered light must provide a substantial contribution to 
their nuclei in both the optical and NIR bands. The lack 
of broad emission line in HIGs is universally associated 
with a substantially decreased luminosity of their nuclei 
in the infrared band and to a different angle of sight to 
the nucleus. This implies that obscuration still plays a 
substantial role. 

The detailed comparison of radio, optical, and infrared 
data for the 3CR sample also shows that the differences 
in the relative contribution in these bands, on a object- 
by-object basis, is dominated by variability and not by 
absorption. Thus, only simultaneous nuclear data would 
allow us to explore intrinsic differences among the sources 
of the sample. The lack of such data effectively ham- 
pers the use of the nuclear Spectral Energy Distribution 
to enlighten how they vary with e.g., source power and 
viewing angle. 
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APPENDIX 

It is interesting to investigate the radio and NIR correlation for FR Is by considering radio data of higher angular 
resolution. To this aim, we take from the literature the fluxes at 5 GHz obtained with VLBI observations. The FR Is 
for which the VLBI radio data are available are listed in Table 5. In Fig. [11] (left panel) we show the correlations in 
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TABLE 5 

VLBI RADIO DATA 



name 


FvLBI 


LvLBZ 


IXGI 


Qf"" Q 1 


OA C\A 

-Z4.U4 




nnc; 
VtUO 


oL- DD 


-Zo. ( 4 


on 1G 

ou.iy 


pnr. 


ov^ oo.l 


OA £^7 


on t^n 


vnn 

Ayy 


54 


OO A A 


Q 1 QO 

ol.oz 


IVUO 


oL> zu4 


OQ Til 


on on 

OU.zU 


pnc; 




9^ 7A 


9Q 97 
zy.z ; 


TQ7 


3C 272.1 


-23.74 


28.67 


G05 


3C 274 


-22.40 


30.02 


G05 


3C 296 


-24.19 


29.83 


G05 


3C 310 


-24.06 


30.66 


G05 


3C 317 


-23.53 


30.80 


VOO 


3C 338 


-23.98 


30.25 


G05 


3C 346 


-23.66 


31.95 


G05 


3C 449 


-24.43 


29.35 


G05 


3C 465 


-23.61 


30.61 


G05 


NGC 6251 


-23.44 


30.61 


J86 



Note. — Table of radio core flux and luminosity from VLBI observations. Col. (1): name; Col. (2): radio core flux in erg s 1 cm 2 
Hz" 1 ; Col. (3 ) : radi o core luminosity in erg s -1 Hz" 1 ; Col. (4 ) : Refe rences:G05: IGiovannini et al.l ||2005|1. X99 : IXu et al.l <1999bf). K05: 
IKovalev et "all (t2005h (VLB A at 15GHz), J97: IJones fc Wehrlel fL997h (VLBA at 1.6GHz). VOO: IVenturi et al.1 020001 ), J86: IJones et all 
111989) . 




the Lvlbi-Lir planes. The linear correlation coefficient is r = 0.90 (P =4.1 xl0~ 6 ). These values as well as the slope 
and the intercept of the derived correlations are essentially identical to those obtained with the VLA data. The reason 
of these similarities can be explained by looking at Fig. [TT] (right panel), which shows the VLA core flux of FR Is 
plotted vs. their VLBI flux. The figure clearly shows that the two measurements lead to very similar flux values. 
The dispersion is extremely small, with the exception of only two outliers (3C 84 and NGC 6251), for which the VLA 
flux is higher than that of VLBI by a factor of ~ 2.5 and ~ 12, respectively. However, since these two objects have a 
relatively high core dominance, they are likely to be observed at a small angle with respect to jet axis and thus they 
are more likely to be affected by flux variability. That might explain the larger discrepancy between the VLA and 
VLBI fluxes for these two sources with respect to the other FR Is. Summarizing, the correlation between the VLBI 
and the IR nuclear luminosity is statistically indistinguishable from that obtained using the VLA cores simply because 
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the core fluxes measured with VLA and VLB1 are only marginally different. 
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